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Abstract. We present a 1+3 covariant discussion of the contribution of gravitational 
waves to the anisotropy of the cosmic microwave background radiation (CMB) in 
an almost-Friedmann-Robertson- Walker (FRW) universe. Our discussion is based in 
the covariant approach to perturbations in cosmology, which provides a physically 
04 ■ transparent and gauge-invariant methodology for CMB physics. Applying this 

approach to linearised gravitational waves, we derive a closed set of covariant equations 
i^^ ' describing the evolution of the shear and the Weyl tensor, and the angular multipoles 

of the CMB intensity, valid for an arbitrary matter description and background spatial 
>0 \ curvature. A significant feature of the present approach is that the normal mode 

^jT ' expansion of the radiation distribution function, which arises naturally here, preserves 

^\ . the simple quadrupolar nature of the anisotropic part of the Thomson scattering source 

o 

& 

in 
c3 



terms, and provides a direct characterisation of the power in the CMB at a given 
O ,■ angular multipole, as in the recently introduced total angular momentum method. We 

provide the integral solution of the multipole equations, and analytic solutions for the 
£h ' shear and the Weyl tensor, for models with arbitrary spatial curvature. Numerical 

results for the CMB power spectrum in open models are also presented, which provide 

an independent verification of the calculations of other groups. 



H ! PACS numbers: 98.80-k, 98.70.Vc 



1. Introduction 



The cosmic microwave background radiation (CMB) plays an essential role in modern 
cosmology. Improved observations of the structure in the CMB should allow one 
to distinguish between rival theories of structure formation, and, in the case of the 
inflationary scenario, to determine the spectrum of primordial scalar (density) and 
tensor (gravitational wave) perturbations. Vector perturbations are expected to leave a 
significant imprint on the CMB only in the presence of active seeds, such as the defect 
models considered in [|TJ. 

The calculation of CMB anisotropics in specific cosmological models is now a large 
industry. For a representative sample of the literature dealing with the imprint of 
scalar perturbations, see e.g. @-|§; for tensor perturbations see e.g. [[J-[]nj. Early 
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calculational schemes were based on linear perturbation theory in the synchronous 



gauge P0[ , with Legendre expansions of the relativistic distribution functions || ||]. 
Some technical improvements in these schemes resulted from the use of gauge-invariant 
variables || [Hj], and, more recently, from improved normal mode expansions of the 
distribution functions |22], p3"] . Adopting gauge- invariant variables frees one from the 
need to keep track of the implications of any residual gauge freedom in the calculation 
(the so called gauge problem; see e.g. |24|, p5fl), while improvements in the normal mode 
expansions have led to more direct characterisations of the anisotropy in vector and 
tensor modes, and a simpler form for the anisotropic Thomson scattering terms than is 
possible with the traditional Legendre expansion. 

Recently, the 1+3 covariant and gauge-invariant approach to perturbations in 
cosmology, developed by Ellis and coworkers [^4), [26J from the covariant hydrodynamics 



of Ehlers P7|] , has been applied to problems in CMB physics p8fl-|3q|. The 1+3 
covariant approach provides a very convenient framework within which to study the 
exact intrinsic dynamics of cosmological models. The approach is formulated at the 
level of variables that are covariantly defined in the real universe, and are (exactly) 
gauge-invariant by construction, thus freeing one from many of the limitations of 
previous methodologies. Furthermore, the approach admits a covariant and gauge- 
invariant linearisation that allows linearised calculations to be performed in a very direct 
manner. In |32j a general 1+3 covariant scheme for the calculation of linearised CMB 



anisotropies was developed (see also [[!3|, |36[). A complete, closed set of equations, that 
fully describe the linearised evolution of inhomogeneity and anisotropy in an almost- 
Friedmann-Robertson- Walker (FRW) universe with arbitrary spatial curvature, were 
presented. These equations were subsequently extended to include non-linear terms 
B3J , where some qualitative implications of non-linear effects were also noted. A 



in 



significant feature of the covariant approach is the unified way in which scalar, vector 
and tensor perturbations are handled; the general perturbation equations are formulated 
without recourse to a decomposition into scalar, vector and tensor modes, although 
such a decomposition can still be made, if required, at a late stage in the calculation. 
In |32| the detailed equations describing scalar modes were derived from the unified 
equations, along with the necessary machinery to calculate the power spectrum of the 
CMB temperature anisotropy from a solution of these equations. It is our purpose here 
to report the extension of the 1+3 covariant analysis to tensor modes. We shall see that, 
in addition to the advantages of gauge-invariance and physical clarity that are evident 
in the scalar calculation, the tensor calculation highlights further advantages of the 
covariant approach over more traditional treatments. Namely, in the present approach 
a normal mode expansion of the radiation distribution function arises naturally that 
has the virtue of maintaining a direct characterisation of the power in the CMB at 
a given multipole £, and of preserving the simple quadrupolar nature of that part of 
the Thomson scattering source terms that arises from the angular dependence of the 
Thomson cross section. Similar advantages are shared by the recently introduced total 



angular momentum method [£2|, [23[ . 
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The paper is arranged as follows. In section |2] we briefly review the 1+3 covariant 
approach to cosmological perturbations in an almost-FRW model, radiative transfer, 
and the decomposition of the CMB sky in angular multipoles. Section [3] begins with a 
discussion of the covariant description of gravitational waves in cosmology, and presents 
linearised wave equations, valid for a general matter stress-energy tensor, for the electric 
and magnetic parts of the Weyl tensor and the shear. These equations extend the 
perfect-fluid equations given in |37|. Following an expansion in tensor harmonics, in 
section § we give analytic solutions for the shear and the Weyl tensor during matter 
and radiation domination, extending the solutions given in [|37| to the case of non- 



flat spatial sections. We also provide the integral solutions to the mode-expanded 
radiation multipole equations in the general case. In section |5| we derive the CMB power 
spectrum in terms of the mode-expanded multipole moments, and present numerical 
solutions in open CDM models, verifying the results in |18 1 . We close in section |6] with 



our conclusions. |Appendix A| contains useful new results on the symmetric trace-free 



representation of the tensor harmonics which we use in sections |] and |5|. 

We employ standard general relativity with a (H ) metric signature. Our 

conventions for the Riemann and Ricci tensors are fixed by 2V[ Q , Vb\u c = R a bc d Ud, 
and R a b = R ac b ', where V a denotes the covariant derivative. Round brackets denote 
symmetrisation on the enclosed indices, square brackets denote antisymmetrisation, and 
angle brackets denote the projected symmetric trace-free (PSTF) part. We use units 
with c = G = 1 throughout, so that the constant in the Einstein equation is n = 87r. 
It should be noted that our convention for the metric signature, which coincides with 
that in our earlier work on CMB anisotropies in the 1+3 covariant approach |3T], |32|, [35 



differs from much of the 1+3 covariant literature. It is straightforward to transform from 
our conventions to those in El: g a b — > —g a b, V a — > V a , and so R a bc d -^ R a bc d - It follows 



that R a b — > -R a fe, but for the Ricci scalar R — > —R. The matter stress-energy tensor 
T a b is unchanged, as is the fundamental velocity (vector) field u a of the cosmological 
model. Raising or lowering an index incurs a sign change when transforming between 
signatures. The transformations of other key geometric objects are discussed in the text 
as they are introduced. 

2. The 1+3 covariant approach to perturbations in cosmology 

The essence of the 1+3 covariant approach is to identify a set of covariant variables 
which conveniently describe the inhomogeneity and anisotropy of the universe in a 
physically transparent (gauge-invariant) manner. Having chosen such variables, exact 
equations governing their dynamics can be found from the full field equations. Since the 
variables are chosen to have a direct physical interpretation, the physical content of their 
dynamical equations is often simple to extract. If required, the exact equations can be 
linearised around a chosen background to yield a simple set of perturbation equations 
that describe the evolution of physically relevant measures of the inhomogeneity and 
anisotropy of the universe in a very direct fashion. 
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The approach employs a 1+3 decomposition of geometric quantities with respect 
to a (fundamental) timelike velocity field u a . We require that, in a general cosmo logical 
model, u a be chosen in a physical manner such that in the FRW limit it correctly reduces 
to the preferred velocity. This restriction is necessary to ensure gauge-invariance of the 
approach. Tensors are decomposed irreducibly with respect to u a into projected tensors 
which are orthogonal to u a on every index. For example, the matter stress-energy tensor 
Tab = T( ab) decomposes as| 

Tab = pu a u b + 2u( a qb) - phab + n ab , (2.1) 

which defines the energy density p = T ab u a u b , the energy flux q a = h b a T bc u c , the isotropic 
pressure p = —h ab T a b/3, and the anisotropic stress n ab = T^ ab ) which is the projected 
symmetric trace-free (PSTF) part of T ab . Here, h ab = g a b — u a u b is the projection tensor 
into the instantaneous rest space of u a (and g ab is the metric of spacetime). In the FRW 
limit q a and n ab vanish since the stress-energy tensor is then restricted to perfect-fluid 
form. 

It is convenient to introduce a projected covariant "derivative" D a which acts on a 
general tensor T b f'" according to 

B a T b J- = h d a h e b ...h c f ...V d T e J-. (2.2) 

It should be noted that D a is a derivation only when acting on projected objects. When 

u a is irrotational (see below), D a reduces to the covariant derivative in the hypersurfaces 

orthogonal to u a . It is also convenient to introduce the derivative along the flow lines 

of u a , denoted by an overdot, so that T b c - = u a V ' a T b .. c '" . With these definitions, the 

covariant derivative of u a decomposes as|]] 

1 
V a u b = u a A b + Gab + -®h ab + u ab , (2.3) 

which defines the acceleration A a = u a , the shear a ab = D( a u b ), the (volume) expansion 
scalar O = D a u a = 3H, with H the local Hubble parameter, and the vorticity 
uj a b = D[ a Wfe]- Introducing a projected, totally antisymmetric tensor e abc = rj abc dU d , 
where ri abc d is the totally antisymmetric tensor on spacetime, we can define the vorticity 
(projected) vector as the dual of u ab \ u> a = e a b c u bc /2. It is also convenient to generalise 
the notion of the 3-dimensional curl to tensors T ab ^, c via 

curlT afe ... c = e de(a D d T 6 ... c) e , (2.4) 

which is PSTF for T ab ^ c PSTF. The vorticity vector can then be written as uo a = 
curlw a /2. The kinematic variables A a , uo a and a a b characterise anisotropy in the local 
expansion and thus vanish identically in an exact FRW universe. The projected gradient 

§ To transform to the conventions of J34], one should make the replacements q a — > —q a , ^ab ~ * ^ab, as 
well as u a — > —u a and h a b — ► —h a b- Of course, p and p are unchanged. 

|| To transform to the metric signature of J34j, make the replacements A a — * —A a , a ao — > —cr a b, 
tOab — > Lu a b, as well as u a — * — u a , V a — > V a and h a b —* —h a b- It follows that for the vorticity vector 
Lu a — ► <-o a - Our Lu a is physically minus that in p4| (rather than just any difference due to the metric 
signature), but coincides physically with that in 
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of the expansion scalar, D a G, describes inhomogeneity in the expansion. Like the 
projected derivative of any covariant object, D a O vanishes in the FRW limit. 

The locally free part of the gravitational field is described by the Weyl tensor C a bcd- 
The ten degrees of freedom in C a bcd can be encoded in two PSTF tensors, E ab and H ab , 
denoted the electric and magnetic parts, respectively:^ 

E ab = u c u d C acbd (2.5) 

H a b = -tacdPbe U & . (2.6) 

The electric part of the Weyl tensor plays an analogous role to the tidal tensor derived 
from the potential in Newtonian gravitation, while H ab is essential for the propagation 
of gravitational waves. Indeed, the non- vanishing of the super-energy flux [35| P a = 



e ac d E bc H b d for any timelike u a forms the basis of Bel's first condition for the existence 
of gravitational radiation in an (exact) solution of the field equations (see e.g. pOfl ). 
Since the Weyl tensor vanishes in the FRW limit, E ab and H a b provide a covariant and 
gauge-invariant characterisation of the perturbation in the free gravitational field. 

To construct a covariant linear perturbation theory about an FRW universe, those 
covariant variables that vanish in an exact FRW universe are regarded as being 0(e) in a 
smallness parameter e. The linearisation procedure then consists of neglecting products 
of 0(e) variables in the dynamical equations. The exact dynamical equations, which 



follow from the field equations, are given in |34], [38|]. Here, we shall only require the 

equations in linearised form around an FRW universe. There are seven propagation 
equations describing evolution along the flow lines of u a : 

- D a q a (2.7) 

±K(p + 3p)+D a Aa (2.8) 

-(p + p )A a + D aP -D b 7r ab (2.9) 

f-curlA Q (2.10) 

- E ab - -nir ab + D {a A b) (2.11) 



p = 


-{P + P, 


e = 


4 
-3**.- 


Qa = 


Ua = 


2 

--Q^a 


0~ab = 


2 

- -QcFab 



Eab = - ®E ab + curliJ ab + -k 



■(p + p)cr ab - D( a g fe) + -hab + -0vr afe 



(2.12) 



H a b = - QHab - cur\E ab - -KCUrl7T afe , (2.13) 

and five constraint equations: 

D a uj a = (2.14) 

DVo6 - curler - ^D,6 - nq b = (2.15) 

D a E ab - k (\eq b + ]-D b p + -D a 7r a6 ) = (2.16) 



,3 J 3 ' 2 

If Replace E ab — ► -E ab and H ab —> -H ab , and also C acbd —> -C ac bd, u a — > u a , and e acd —> e acd to 
transform to the conventions of [Ml . 
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D a H ab - -K[2(p + p)u b + curlg b ] = (2.17) 

H ab - curlcr a6 + D( a a; b ) = 0, (2.18) 

which serve to constrain the initial data. It is straightforward to verify that the 
constraint equations are preserved by the propagation equations. 

We describe the CMB radiation field via the specific intensity I(E, e a ) which gives 
the intensity per unit solid angle along the (projected) direction e a at photon energy E. 
The energy and propagation direction relative to u a are related to the four-momentum 
of the photon by 

p a = E{e a + u a ). (2.19) 

The photon distribution function I(E,e a )/E 3 is frame-independent for given p a . It is 
convenient to work with the energy-integrated intensity, 



I(e a ) = / I(E } e a )dE, (2.20) 

Jo 

since this determines the bolometric temperature along a given direction e a . Denoting 

the fractional temperature difference from the all-sky average at some point by (5r(e a ), 

we have 

[l + Me a )] 4 =(^//(e a )dfi) '/(e-), (2.21) 

where the integration is over solid angles. 

The anisotropy of the CMB intensity is conveniently described by the PSTF 
multipole moments i^j^in the covariant angular decomposition of I(e a ): 

oo 

I(e a ) = J2Aj x I Ai e A *, (2.22) 

e=o 

where e Ae is a convenient condensed notation for e ai . . . e ae , and we have defined 
_ 4tt(-2)'(1!)* 

Ae = (2£ + iy. ■ (2 - 23) 

Including the factor of A^ in (|2.22|) has the virtue of ensuring that the three lowest 
multipoles /, I a and I ab are simply the radiation energy density, energy flux, and 
anisotropic stress, respectively: 

p (7) = j 5 g W = J ft) l$ = I A . (2.24) 

Propagation equations for the radiation multipoles follow from a multipole expansion 
of the Boltzmann equation |32|, |34|. In linearised form, we have 

4 h £ 4 8 

= - n e a T (l Ai - 8° e I - is}lv® - ^|/ 0ia2 ) , (2.25) 

with Ia,, = for £ < 0. The projected vector v^ which appears in the scattering 
term on the right of ( |2.25| ) is the relative velocity of the baryon frame, so that 

+ We use a lumped index notation, Ai = a± . . . a^. 
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v$ = u a + v^ + 0(e 2 ) where u^ is the baryon four- velocity. The electron number 
density in the baryon frame is denoted n e , and <tt is the Thomson cross section. We 
have ignored the small effect of polarization in ( |2.25| ). In \4 i\ we give a full 1+3 covariant 



treatment of CMB polarization. There it is shown that the only effect of polarization 
on the intensity multipole equations is to replace I ai a 2 by Ia X a 2 + 6£ a ia 2 hi the final term 
on the right-hand side of ( |2.25| ), where S aia2 is a rank-2 PSTF tensor which describes 
the quadrupole moment of the electric polarization; see [^| for equivalent results in 
the total angular momentum approach. It is clear from the multipole equations that 
the complications due to the anisotropy of Thomson scattering affect only the £ = 2 
multipole equation. (This statement is also true of the further complications due the 
polarization dependence of Thomson scattering.) This simplicity, which is preserved by 
the normal-mode expansion described in section [| ensures that the current methodology 
shares many of the advantages of the total angular momentum method over traditional 
Legendre techniques. 

It follows from equation (|2.21|) that in linear theory, <5r(e a ) has the multipole 



decomposition [28 



MO^E^'V- (2-26) 

1 £=l 

so that in a statistically isotropic ensemble the power spectrum of temperature 
anisotropies, Ci, where 

00 IP + 1 
= 52^-—C e Pe(cos6), (2.27) 

l=\ 47r 
be determined from the covariance of the radiation 





(6 T (e a )S T (e' a )) 


with cos 9 


= —e a e' a , can 


multipoles 


33: 




(j)\w B *) 


Here fcgjj s 


= h,* . . -h n e \, anc 

(ai at) ' 



AeCeSet'h^y (2.28) 

and large angle brackets denote an ensemble average. 

3. Gravitational waves in an almost-FRW universe 

The 1+3 covariant description of gravitational waves in a cosmological context has 
been considered by Hawking [[42], and more recently by Dunsby, Bassett, & Ellis |37|| . 



Linearised gravitational waves are described by the transverse degrees of freedom in the 
electric and magnetic parts of the Weyl tensor. Similarly, the shear and anisotropic 
stress are transverse so that 

B a E ab = 0, D a H ab = 0, DV o6 = 0, D a ir ab = (3.1) 

at 0(e). Furthermore, the vorticity and all projected vectors that are generally 0(e) in 
an almost-FRW universe vanish at linear order. In particular, in a pure tensor mode 
the individual matter components all have the same four-velocity which defines the 
fundamental velocity u a . 
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The covariant equations governing the propagation of linearised gravitational 
waves can be obtained from the equations of the previous section. The constraint 
equations fl2.14J - p7TB ) reduce to the single equation 

H ab = curla a6 . (3.2) 

which determines the magnetic part of the Weyl tensor from the shear. The propagation 
equations fl2.7H2.13D reduce to 

p = _( P + P )e (3.3) 

(3.4) 
(3.5) 
(3.6) 
(3.7) 



e = 


~e 2 -^(p+3p) 


O'ab = 


2 1 

— 77 ® a ab ~ E a b — y K7V ab 


E a b = 


- QE ab + curl H ab + -k 


H a b = 


- QH ab - cur\E ab - -k( 



■{p + p)o a b + T^ab + ^@7r a6 



Making use of the linearised identities p9| , valid for an arbitrary first-order PSTF tensor 

^abi 



u b V b ( curl Sab) 



curl5 a6 - -6 curl S a 



D a S, 



ab 



:Q-D a S a 



-curl(D a 5 a 



u c \7 c (D a S ab ) 

D a (cur\ Sab) 

it is straightforward to verify that the constraints in equations ( |3.1|) and 
consistent with the linearised propagation equations ( |3.3| -|3.7|). 

Inhomogeneous wave equations for the shear and the magnetic part of the Weyl 
tensor follow from differentiating equations (|3.5|) and fl3.7|) along the flow lines of u a : 



(3.8) 

(3.9) 

(3.10) 
2|) are 



(Tab + D 2 0"a6 + "©^ab + 

H ab + D 2 H ab + 7 -QHab + 2 



1/ x K 

2(4-37)^-^ 



[2-^Kp- 



3K 



O'ab = —K 
H a b 



^ab + -@7Tab 



— K 



curl7T afe + -6curl7r afe 

o 



where D 2 = D a D a , and we have used the linearised identity 



curl ( curl S. 



,ib) — D S a b + ~^-S a b 



:D(aD c 5'b} c . 



(3.11) 
(3.12) 

(3.13) 



Here, 7 is defined in terms of the total energy density and pressure: p = (7 — l)p, S 
is the scale factor of the universe [defined in a general almost-FRW model by S = HS, 
and D a S = 0(e)], and 6K/S 2 is the zero-order intrinsic curvature scalar on comoving 
hypersurfaces. 

The solutions to the homogeneous wave equations in a K = universe with constant 
7 were given in [J37H - Solutions for K 7^ are given in section |] in the radiation 
and matter eras. Solutions to the inhomogeneous equations can be obtained from 
those solutions by Green's method, but an accurate treatment of the matter-radiation 
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transition requires a numerical integration. The damping due to the anisotropic 
stress terms in the wave equations was discussed by Hawking for the simple viscous 
model ii a b = Xcr a b, where A is the viscosity f£2fl . During tight-coupling, the radiative 
contribution to 7T ab for a general (long wavelength) perturbation is 



l ab 



16 

27 



It c (a a b + D( a u 



b) > 



(3.14) 



where t c = (ricOr) -1 is the photon mean free time. Including polarization changes 
the factor of 16/27 to 32/45 in this expression for I ab , since the electric polarization 
quadrupole £ ab ~ / a b/4 during tight-coupling. The tensor a a b + D( a ?A is the shear in 
the baryon frame, so the effective (polarization-corrected) viscosity is A = 32Jt c /45. 
Radiative viscous effects, as well as those due to other (collisionless) species such as 
neutrinos, lead to a small increase in the damping of gravitational waves over that due 
to redshift effects alone. However, unlike the case of scalar perturbations, the photon 
(and neutrino) quadrupole does not leave a significant imprint on the tensor CMB 



temperature power spectrum at high multipoles |19|| , since on the scales where damping 
due to anisotropic stresses could be significant, the tensor power spectrum is already 
insignificant due to sub-horizon redshifting of the gravity waves. 

The wave equation for the electric part of the Weyl tensor is not closed (even for 
K a b — 0) because of the presence of terms involving the shear: 



E ab + D 2 E ab + -QE ab + 2 
1 



3K 

(2- 7 )kp- — 



Eab + -I^Pl 

6 



- -KD 7T ab + 



-np 



3K 



K1X ab 



6 



-Qktt, 



ab 



3- + (2 - 3 7 )6 

7 



-KKab- 



(Tab 



(3.15) 



To eliminate the shear requires a further differentiation j|2| (the resulting equation for 
T^ab = was given in |37|). Although one cannot write down a closed wave equation 
for E ab , the effective order of the equation is two since E ab can be determined from the 
shear via equation ( |3.5| ), and the shear is determined by a wave equation. 



3.1. Mode expansion in tensor harmonics 

In linear perturbation theory it is convenient to expand the 0(e) variables in harmonic 
modes, since this converts the constraint equations into algebraic relations and the 
propagation equations into ordinary differential equations along the flow lines. For 
tensor perturbations, the appropriate harmonic functions can be derived from the 
tensor harmonics Q ab . The zero-order properties of the tensor harmonics are discussed 
in [Appendix A], where their explicit PSTF representation is derived for general K. 



The electric and magnetic parts of the Weyl tensor, the shear, and the anisotropic 
stress can be expanded directly in the electric and magnetic parity tensor harmonics 
(see Appendix AT ): 

(3.16) 
(3.17) 



E ab =S- 2 J2k 2 (E k Q a k b ) +E k Q 



H, 



)(*) 



ab i 



(^ 



„/, = S ^k (H k Q a l' + H k Q ab 
k 
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(3.18) 

(3.19) 



where the symbolic J2k denotes a sum over the harmonic modes. For closed spatial 
sections (K > 0), the spectrum is discrete with v = \{k 2 + 3K)/|K|] 1//2 an integer 
> 3. We use an overbar to distinguish the magnetic parity harmonics from the electric 
parity, and a subscript/superscript k to label (implicitly) the distinct eigenfunctions with 
degenerate eigenvalue k. For the representation in |Appcndix A , the label k represents 
the collection v, £, m, where £ and m describe the orbital angular momentum of the tensor 
harmonic. Functions of the eigenvalue alone will be denoted with k as an argument, 
for example, A(k), or, equivalently, A{y). The mode coefficients, such as E k have 0(e 2 ) 
projected gradients by construction. 

For the specific representation of the tensor harmonics constructed in |Appendix A| , 
the electric and magnetic parity tensor harmonics are related by a curl: 



curlQ 



ab 



k / 3K\ 1/2 

s( 1 + ^) 



6 {k) 



(k) 



with an equivalent result for curl Q ab 
in a k : 



From equation 



Hh 



3if\'/ 2 



' w) 



0~k- 



(3.20) 
we find that H k is algebraic 

(3.21) 



Note how the presence of curl terms leads to a coupling of the different polarization 
states. The constraint equation ( |3.2| ) [or, equivalently, fl3.21p ] allows one to eliminate 
H a b from the discussion, leaving coupled first-order equations for the electric part of the 
Weyl tensor and the shear: 



k 2 /* l n „ \ k 

p [Bu + - 3 e Et ) - - 
s {"" + 3 eak ) + T* Ek 



'P 3K_ 

i 



i 



i 



- ;t;p\ o- fc = -g(37 - l)«p©7T fe -I -h-f)Ki, 



--^P^k- 



(3.22) 



(3.23) 



Combining these equations gives the harmonic expansion of the wave equation for the 
shear [equation (|3.11|) ]: 

a k + Qa k + 7^7 + ^7 - -V6"f-2)Kp a k 



S 
KP 1 



-(37 
3 V ' 



2)On k - 7r k 



For the radiation multipoles, Ia v we expand in PSTF tensors Q 
overbarred counterparts) derived from the tensor harmonics via 

' ' x D/ D {k) 



(k) 

At 



Q 



(k) 



The factor of (S/k) 1 2 ensures that Q A = 0(e). For the £-th multipole, 



k 



1 k ^Ati 



(3.24) 
(and their 

(3.25) 
(3.26) 
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where, for later convenience, we have defined 



Pi = n *«> 



(3.27) 



n=2 



K t = [i- (r - 3)K/k 2 ] 



211/2 



> 2. (3.28) 

I, — K-2^ 1 ll ■ The Q;J satisfy a number 
of differential identities (which are independent of any specific representation for the 
tensor harmonics) , the most useful of which are 



It follows that the photon anisotropic stress ir^ 1 



B ae Q 



At 



(£ 2 -4) k 



d 2 qS 



£(2£ - 1) 

s 2 



s 



l-(f 






Q { t 



-2)(* + 3) 



K 



Q 



(k) 

At 



for C. > 2, and the recursion relation 

k. 



Q 



(k) 



g^iatQAt-i) 



(3.29) 
(3.30) 

(3.31) 



>(*) 



for £ > 3. In a closed universe, the Q^ vanish for £ > z/, so that only modes with v > £ 
contribute to Ia v The distinct processes of decomposing J(e a ) into angular multipoles 
Ia v followed by an expansion of the multipoles in derivatives of the tensor harmonics, 
combine to give a normal mode expansion of the radiation distribution function which 
may be regarded as the extension of the Legendre tensor approach, first introduced by 
Wilson ||, to non-scalar perturbations. In our opinion, it is one of the virtues of the 1+3 
covariant approach that the physically significant decomposition into angular multipoles 
I At, is kept distinct from the mathematically convenient harmonic mode decomposition. 
Maintaining this distinction allows one to describe the evolution of the anisotropy in 
a general cosmological model directly, without recourse to a harmonic splitting of the 
perturbations. Moreover, the present treatment is manifestly gauge- invariant, although 
we do not labour this point since gauge-invariance is less of an issue for tensor modes 
than for scalars. 

Inserting the mode expansion, equation (|3.26| ), into the propagation equation (|2.25| ) 
for the multipoles Ia ( , we find the mode-expanded multipole equations: 



/?> + 



(l + 3)(l-l) 

(£ + l)(2£+l 



-Kl+l+k 



(*+l) 



(2£ + l 



-Ktl 



(£-1) 



8 k 



r2 

OpK 2 ak 



n e a T I, 



r(0 



2r(2) 



To 6 '' 



l x k 



(3.32) 



with similar equations for the barred variables. Massless neutrinos satisfy the same 
multipole equations but with the Thomson scattering terms omitted. If polarization 



r(2) 



r(2) 



(2) 



is included, I k should be replaced by I\. ' + Q£j, in the last term on the right-hand 

/r)\ 

side of (|3.32|), where E\ is the coefficient in the harmonic expansion of the electric 



polarization quadrupole; see |^TJ for further details. Equations ( |3.24j ) and ( f3.32[ ), plus 
the massless neutrino equations, form a closed set of equations that can be solved to give 
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a complete description of the tensor contribution to the CMB anisotropyQ. Analytic 
solutions for the shear and electric part of the Weyl tensor are given in the next section, 
assuming that anisotropic stress effects are negligible, as well as the integral solution for 
the radiation multipoles. 

We end this section by noting that in kinetic theory, where the radiation multipoles 
may be non-vanishing for all £, the familiar decomposition of the perturbations into 
scalar, vector and tensor modes is incomplete. A complete set of eigenf unctions for 
expanding a rank-£ PSTF tensor must include the equivalents of the scalar (rank-0), 
vector (rank-1), and tensor (rank-2) harmonics up to rank-£. Since the perturbations 
in the matter stress-energy tensor, the kinematic variables, and the Weyl tensor are 
not affected by the higher-rank perturbations in linear theory (where perturbations of 
different rank evolve independently) the higher-rank perturbations in the distribution 
function cannot arise spontaneously. Any power that was initially in higher-rank 
perturbations of the photon distribution would thus be damped away due to Thomson 
scattering, leaving a negligible effect on the present day CMB anisotropy. 

4. Analytic solutions for E a b, a a b and I^ t 

During tight-coupling, the photon anisotropic stress is suppressed by Thomson 
scattering, so the only significant contribution to TT a b is from the neutrinos. Assuming 
that this contribution is negligible, it is straightforward to find the analytic solutions 
for the Weyl tensor and the shear during radiation or matter domination. Anisotropic 
stresses can be included analytically by constructing the Green's functions from the 
solutions of the homogeneous equations given here. 

It is convenient to use x = y/\K\rj as the independent variable, where rj is 
conformal time. During radiation domination in an open universe, HS = y/\K\ cothx, 
so that (|3.24j) becomes 



d <*k , n „ da fc , / 2 , 1 2 



. , + 2 cothr— - + [i/ + 1 - — — a k = 0. (4.1) 

dx z ax V smh x) 

The solution of this equation which is regular as rj — > is 

cosz/x sinz/xcothx 

o- k oc — — , 4.2 

smhx v smhx 

and the irregular solution is 

sin vx cos vx cothx 

(Jk oc — — + — . 4.3 

smhx v srnhx 

These solutions for the shear determine the electric part of the Weyl tensor from 



equation (|3.23| ). The regular solution evaluates to 

1 



E k oc 



\Jv 2 + 3 sinhx 



smra , 2 

[y + 1 — coth x) + cosz/xcothx 



(4.4) 



v 

* When polarization is included, additional multipole hierarchies must be included for the electric and 
magnetic multipoles of the polarization to close the equations JO] . 
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and the irregular solution to 

E k oc 



\/v 2 + 3 sinhx 



cos vx , 2 « i 2 \ • n 
(v + 1 — cotn x) — sin vx cothx 



(4.5) 



v 

For kt] <C 1 (which implies vx <ti 1 and x <C 1), the mode is well outside the 
horizon since kr] <C 1 implies that HS/k m 2/(ki]). When this condition is met the 
electric part of the Weyl tensor is approximately constant, with equation (fOJ) giving 
E k oc 2(v 2 + l)/[3y/(v 2 + 3)]. The solutions for a closed universe can be obtained from 
those in an open universe with the substitutions x — > ix and v — > —iv. 

During matter domination in an open universe, HS = y/\K\ coth(x/2), so 
that (|3.24D becomes 

^ +2 coth( I /2)^ + (^ + l-- i? L_).^0. ,4.6, 

The regular solution of this equation is 

I fsini/T 1 

o k oc 2/ ; - <^ [Av 2 + 1-3 coth 2 (a;/2)] + 3 cos vx coth(z/2) \ , (4.7) 

sinh (x/2) I 2v J 

and the irregular solution is 

I ( cos &% 1 

a k oc 2/ - - [Av 2 + 1 - 3coth 2 (x/2)] - 3 sin vx coth(x/2) \ . (4.8) 

sinh (x/2) I 2v J 

Using equation ( |3.23| ), the regular solution for the electric part of the Weyl tensor during 
matter domination is 

E k ex X {^^coth(x/2)[2^ + 1 - coth 2 (x/2)] 

V^ +3smh (x/2) ( 2v 

- cosz/x[2z/ 2 + 2-3 coth 2 (x/2)] I, (4.9) 



and the irregular solution is 

1 j 3 cos vx 

V^ 2 + 3 sinh 2 (x/2) \ 2v 

+ sinz/x[2z/ 2 + 2-3coth 2 (x/2)] I. (4.10) 



E k oc ; = == , i 2 — ^ — ^^ coth(x/2)[2i/ 2 + 1 - coth 2 (x/2)] 



On superhorizon scales, E k oc 2(i/ 2 + l)(Av 2 + l)/[5^/(^ 2 + 3)] is approximately constant 
in the regular mode. In both the matter and radiation eras, a k = 0(x) in the regular 
modes on superhorizon scales. It follows from ( |3.21|) that on superhorizon scales the 
Weyl tensor is dominated by the contribution from its electric part. The solutions given 
here generalise the K = solutions derived in [37]; in the limit K — > we reproduce 
their results. 

We can provide a formal solution to equation ( p.32|) for the intensity multipoles 
by noting that the homogeneous equations (obtained by setting n e = and a k = 0) 
are solved by the functions £(£ — l)$^(x)/sinh (x) in an open universe, where x = 
y/\K\(rjR — if) with rjR the conformal time at our current position R. Here, $£(x) are 
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the ultra-spherical Bessel functions ||; see |Appendix A| also. Defining the optical depth 
back to cosmic time t, with t — tn at R: 



f 
r(t) = / n e a T dt , 



the integral solution for I k 



m an open universe is 



r(Z) 



4£{e 



tR 



dt'e~ 



k 



5 afc + -^ n ^^k 



(7) 



M 



,.r 



sinh x 



(4.11) 



(4.12) 



fc [(z/ 2 + l)(z/ 2 + 3)] 1 /2 

For a closed universe, the hyperbolic functions should be replaced by their trigonometric 
counterparts and z/ 2 + n by z/ 2 — n, where n is an integer. Equation ( [4.12 ) can be shown 
to be equivalent to the tensor integral solution derived in pi with the total angular 
momentum method. The shear term in ( [4.12| ) arises from the cumulative redshifting 
effect of shear along the line of sight. The contribution to the anisotropy at R from 
the redshift incurred for a given increment along the null geodesic is weighted by the 
probability e _r that a photon subsequently does not scatter before arriving at R. The 
term involving the radiation anisotropic stress (intensity quadrupole) in equation (|4.12j ) 
describes the cumulative effect of the anisotropic scattering into the beam along the 
line of sight. (Inclusion of polarization gives an additional contribution to this term 
so that TTj? should be replaced by tij? + Qk^S^ . This term has only a very small 
effect on the intensity multipoles over the range of £ where the magnitudes of the Ia e 
are appreciable.) The geometric factor $£(x)/sinh 2 x arises from the projection of 
the contraction Q ab e e e b of the tensor harmonic with the projected photon direction at 
comoving distance x/\/\K\ back along the line of sight. This follows quite generally 
from noting that, for any representation of the tensor harmonics, Q ab e a e b at x down 
the line of sight is given in terms of the Q A e Al at the observation point by 



Qi?e a e b | 



47T 



[(^ + l)(^ + 3)]i/: 



E 



*(* -1) *?(*) 
Aifle sinh x 



VA, e l-R- 



(4.13) 



This result, which can be regarded as the PSTF extension of the Rayleigh plane wave 
result to eigentensors of the Laplacian on a manifold of constant curvature, can be 
verified by comparing the derivatives of Q ab e a e b at x = with those computed from 



( fc )„A» 



the following recursion relation for the Q A e 



( k ) n A t 



e b D b (Q^e 



k 
S 



Q 



(k) 



i-l+l 



(f 



H+\ 



(4£ 2 - 1 



4 ) K 2 n (k) 



M- 



(4.14) 



In a closed universe the sum over 
2<l<v-l. 



in equation ( |4.13| ) is restricted to the range 



5. The CMB power spectrum 



(*) 



>(*) 



To determine the CMB power spectrum from the I\ we require the Q A at the 
observation point R. For this purpose, it is convenient to consider the variation of 
Q ab e e e b down the line of sight for the chosen representation of the harmonics, and to 
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compare this result with ( |4.13| ) which is independent of the specific representation. For 
the representation of the tensor harmonics given in |Appendix A| , we find for the electric 
parity harmonics 



Q%e a e b 



in an open universe, and Q 



2(i/ 2 + 1) 



(£ + 2)! 
(1-2)! 



1/2 



sinh 2 x 



Q-Afi' 



\R: 



(5.1) 



( k )„a„b 
ab 



e e 



for the magnetic parity harmonics. Here, Qa„ is 



a rank-£ PSTF tensor encoding the 2£+l degrees of freedom in the tensor harmonic with 
£ units of orbital angular momentum; see the discussion in [Appendix A.3| . Comparing 
equations (|5.1|) and fl4.13|) , and using the orthogonality of the e^ Al \ we find 

1/2 r(£ + 2)(£+iy 1/2 



Q 



(*)i 

A„\R 



1 

47T 



(v 2 + 3) 



2z/ 2 



-I) 



AiPiQaA 



(5.2) 



for the electric parity tensor harmonics, and Q\,\r = for magnetic parity. Although 
the magnetic parity harmonics do not contribute to the radiation anisotropy at R, they 
do contribute on other integral curves of u a . The contribution of the electric parity 
harmonics undergoes a compensating change to ensure statistical homogeneity. 

Following standard practice we assume that on the scales of cosmological interest 
today there is negligible power in the irregular modes given in section [| during the 
radiation era. This assumption is justified by the highly squeezed quantum state 
predicted at the end of inflation [[|3|, [44J. Since Ek is constant outside the horizon 
for the regular solution, it is convenient to parameterise the radiation multipoles in 
terms of this constant: 

(^ 2 + 3), 



M 



, !-T ie) (v) lim E k , 



(5.3) 



in an open universe. T^\v) is the radiation transfer function, and the ^-dependent 
factor is chosen to conform with treatments based on the metric perturbation variable 
(e.g. |23|). The primordial tensor power spectrum V[y) is defined in terms of the 
covariance matrix for the Ek in the limit krj — ► 0: 

1 {y 2 + 1) 2 V{y) 



(E k Ek') = (E k Ek>) oc 



-Ofcjfc' 



(5.4) 



\K\ 3 / 2 (is 2 + 3) 2 is(v 2 + 1) 
where 8kk' is the symbolic delta function appropriate to the sum over tensor harmonics; 
see |Appendix A| . In the K — > limit, V{y) —>■ k nT for a tensor spectral index n^. Making 
use of equations ( |2.28| ) and ( |5.2| ), we find that the CMB power spectrum evaluates to 

1 (i + 2)(£ + l) r°° udv (^ 2 + 3), 



C, 



-V(u)[T«\u)f 



(5.5) 



16 £{£-l) Jo (z/ 2 + l) 2z/ 2 
in an open universe. For closed models the integral over v should be replaced by a 
discrete sum over integral v > £ + 1. Note that the CMB power spectrum at multipole 
£ depends only on the transfer functions at the same multipole; a feature not shared by 
schemes based on simple Legendre expansions of the distribution functions (e.g. [14]). 
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With these conventions, the contribution to the magnitude (squared) of the electric part 
of the Weyl tensor from superhorizon modes is 

Jk) dhTk {EabE > = Y u^ + 3) V{v) - (5 ' 6) 

For v — ► in open models, the limit V{y) — > is required to avoid an infrared divergence 
in the Cg and the electric part of the Weyl tensor. (In bubble models, V[y) — > v for 

i/<i no 

The Ci due to gravitational waves in the CDM model are plotted in Figure [l] for 
three values of the current-day matter density: Qq = 0.1, 0.4 and 1. The Hubble constant 
h = 0.7 (Hq = lOOh kms _1 Mpc _1 ), the baryon fraction Q^h 2 = 0.0125, and there is no 
reionisation or cosmological constant. The primordial power spectrum is the minimal 
scale-invariant spectrum V[y) = tanh(7rz//2) fl45 ]. For the flat model (Oq = 1), this 



V[y) —* 1 corresponding to n^ = 0. The power spectra were calculated by modifying 
CMBFAST H to solve the 1+3 covariant equations given in this paper. The results 
are in excellent agreement with those given in [[OJ which employed the total angular 
momentum method. Notable features in Figure [I], which are discussed further in |18|| , 



include the angular scaling brought about by changes in Q (variations in the angular 
diameter distance), the curvature cutoff in the large scale CMB power spectrum in low 
density universes (the curve for Q = 0.1), and the enhancement in power on large 
scales, compared to the flat case, for higher values of the density parameter (the curve 
for Qq = 0.4), due to the growth of the shear in the longest wavelength modes once they 
enter the horizon. 

6. Conclusion 

We have developed the 1+3 covariant results necessary to calculate the temperature 
anisotropy of the CMB due to gravitational waves in general almost-FRW models. 
The covariant and gauge-invariant approach is well suited to describing gravitational 
waves in a cosmological context, since the emphasis is placed on the physically relevant 
Weyl tensor and shear, rather than the perturbation in the metric tensor. In the 
covariant calculation the intensity anisotropy is decomposed into irreducible form before 
linearisation around an FRW background, so the approach is also well suited to studying 
non- linear radiative effects (see f3"l|). For the linear calculation, the expansion of the 



radiation intensity in PSTF tensor-valued multipoles ensures that the subsequent mode 
expansions of the radiation variables for scalar, vector and tensor modes naturally 
generalise the harmonic expansions of the rank-0, 1 and 2 tensors which describe 
inhomogeneity and anisotropy in the geometry. The same is true for the description 



of the polarization of the CMB, as we show elsewhere [^, |41[]. The technique should 
also allow a straightforward derivation of the multipole hierarchies for massive particles 
in difficult cases, such as non-flat models with vector or tensor perturbations. 
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Figure 1. Intensity multipoles from gravitational waves in CDM models with h — 0.7, 
Q,\,h 2 = 0.0125, and no reionisation or cosmological constant. The Ci are shown for two 
open models, Qq =0.1 and 0.4, and the flat model. The primordial power spectrum is 
V(y) = tanh(7ri//2). The curves are normalised at Ciq. 
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Appendix A. The tensor harmonics 

The tensor harmonics Q ab are zero-order, second-rank PSTF tensors satisfying 



Y) 2 {k) 






(A.1) 



where the equality holds at zero-order. They are constructed to satisfy the zero-order 
relations 



D a W = 0> 0$ = 0, (A.2) 

which are consistent with ( |A.1| ) at zero-order. With the definition v = [(k 2 +3K)/\K\] 1 ^ 2 , 
the spectrum of regular, normalisable solutions for open and flat background models has 
v > 0, while in closed models v takes integral values > 3. Supercurvature modes in open 
models (see e.g. [f46| ) correspond to > v 2 > —3. For v 2 ^> 1 the mode is well inside 
the curvature radius and k is an effective comoving wavenumber. We label the tensor 
harmonics with a lumped superscript (k) which should be understood to distinguish 



degenerate solutions of ( |A.1|) . 
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In the text we make use of a specific PSTF representation of the tensor harmonics. 
For this representation it is convenient to "coordinatise" the spacetime as follows. 
Starting at some arbitrary point R (which can conveniently be taken to be our current 
location), we construct spatial geodesies with tangent vectors e a (x), the generators of 
which at R are the set of projected unit vectors {e a |#} at R: 

e a D a e fe = 0, u a e a = 0, — = -Se a , (A.3) 

dx 

with x > (the equality holding at R). In the absence of vorticity the curves x a (x) lie 

in the flow-orthogonal hypersurface through R. A projected vector field e a and a scalar 

field x can then be constructed in the neighbourhood of R from their restriction on the 

hypersurface defined by {e a \R,x} y i a 

e^ = 0, X = 0. (A.4) 



The fields e a and x satisfy equations (|A.3|) at zero-order in an almost-FRW universe. 



To the same order, the restriction of e a to the past lightcone through R is the projected 
direction of a photon that passes through R, and x is the zero-order conformal lookback 
time. 

Appendix A.l. Electric parity harmonics 

For the electric parity tensor harmonics we start with the most general electric parity, 



rank-two PSTF tensor [|7| which is constant along the flow lines of u a : 



Qab = Ti(x) (e a e b Q Ce e c * + ~H ab Qc e e Ct 

+T 2 {x)e {a 7i%Qc t e c ^ + T 3 (x)[Q abCe _ 2 e c ^] TT , £ > 2. (A.5) 

Here, Qa € is a rank-i PSTF tensor field satisfying the zero-order equations 

e b B b Q Ae = 0, Q Al = 0, (A.6) 

which determine the tensor field Q Al from its value at the point R, and Ti(x), i = 1,2, 3, 
are scalar functions of x = y/\K\x- The (screen) projection tensor TC ab is 

Hab = Kb + e a e b , (A. 7) 

and the notation [yl a fe] TT denotes the transverse traceless part of the second-rank tensor 
A ab : 

[A ab ] TT = H c a 'H^A ClC2 - l -H ab n c ^A ClC2 . (A.8) 

Demanding that D a Q ab = determines T 2 and T 3 in terms of T\. In an open universe 
we find 

™ = ( £+i)L^ |3ii ' iAr ' wi - (A - 9) 
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Equation ( |A.l ) gives three further equations, which are satisfied by virtue of 
equations ( |A.9| ) and ( |A.10| ) provided that 



d 2 d 

-— rTi(x) + 6cothx-— TAx) + 
dx z dx 



„* + 9- <^- 6 > 



sinh x 



The regular, normalisable solution of this equation is 
T 1 (x) = N(u)- 



$t(x) 



Ti(x) = 0. (A. 11) 



(A.12) 



sinh 2 a; ' 

with v > (subcurvature). The ^(x) are the ultra-spherical Bessel functions defined 
recursively from 
d 



f x *>w 



'cothx^(x) - \v 2 + {£ + l) 2 ] 1/2 $?+i(aO) 



where 



_ smux 
v sinhx 



We choose the normalisation constant N(u) so that 

J dn e a lR dxsmh 2 xQi% 2 Q^ a ^ = ^S w AeQ Al Q' Al u- 2 S(u -v') 
which gives 



N(u) 



1 



„J2(v* + 1) 



(£ + 2)! 
(£-2)\ 



1/2 



(A.13) 
(A.14) 
(A.15) 

(A.16) 



For closed universes, the hyperbolic functions should be replaced by their trigonometric 
counterparts in equations ( |A.9| - |A7L6| ), v 2 + n should be replaced by v 2 — n with n 
an integer, and 5{y — u') should be replaced by 8 UV '. For closed models, the regular, 
normalisable modes have v an integer > 3, restricted to v > £. 



Appendix A. 2. Magnetic parity harmonics 

For the magnetic parity tensor harmonics, which we denote with an overbar, we start 



with the general magnetic parity, rank-two PSTF tensor [47]: 

Q<£ = f l {x)[e cl e c ^-\ a Q h) c i _ 1 e c ^\ TT + f 2 {x)e d e {a e b) d ^Q Ct e c '-\ (A.17) 
with £ > 2. Equations ( |OD and (TAT2]) give 



TAx) 



(£ + 2) smlTo; dx 



[sinh 3 xT 2 (a;)], 



and 



d 2 - d - 

-— T 2 (x) + 4cothx— T 2 {x) + 
dor dx 



z/ 2 + 4 



[£ 2 + £ - 2) 
sinh 2 x 



The regular, normalisable solution of ( |A.19| ) is 



sinhx 



(A.18) 



T 2 (x) = 0. (A. 19) 



(A.20) 
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where v > 0. With the same normalisation as the electric parity harmonics, we find 



N(p) 



e+i: 



iVa r 



l/» + 1 



(£ + 2)! 
(£-2)! 



■nl/2 



(A.21) 



so that N(v)/N(y) = 2i//{£ + 1). In a closed universe we make the same replacements 
as for the electric parity harmonics. 

As a consequence of the first-order identity 

D 2 f cm\S ab ) = curl (D 2 S ah ), (A.22) 



curl(D S a b), 

where S a b is an 0(e) PSTF tensor, it follows that if Q ( ab solves equation (|A.1|) then so 
does cur\Q ab . Furthermore, by virtue of equations (|3.8| ) and (|3.10|) , cur\Q ab will also 
satisfy (|A.2|) if Q ( ab does. Since the curl operation is parity reversing, we expect the 
electric and magnetic parity tensor harmonics for given £ to be related by a curl. With 
the normalisation conventions adopted here, it is straightforward to show that 

k 
S 



curl Q 



(*) 



3K\ 1/2 -, 



t;(1— S Q<$, 



/.' / 3^\ 1/2 



curlQi? = ^l + ^j Qi\ 



which are consistent with the first-order identity ( |3.13| ) 



(A.23) 

(A.24) 



Appendix A. 3. Sum over harmonic modes 

A normalisable, first-order, rank-2 PSTF tensor field, such as the electric part of the 
Weyl tensor, can now be expanded in terms of the electric and magnetic parity tensor 
harmonics. Since the Qa £ have only 2£ + 1 degrees of freedom, it is convenient to 
introduce a set of 2£ + 1 orthogonal basis PSTF tensors {Qj\T }, with m = —£, ...,£. 
We choose the basis tensors so that 



q(*™) Q(tm')At _ a-U 



which implies 



Y,Qt ] Q {em)Be 



-lu(Bt) 



Aj'h 



e ,0 (A e )- 



(A.25) 



(A.26) 



The tensor harmonic expansion of E a b can then be written as 



E ab = / \K\ 3 / 2 u 2 du J2 E (KimQZ + E vtm Q$), (A.27) 

Jo e=2 m=-e 

in an open universe. [Equation (|A.27|) can easily be generalised to include supercurvature 

modes, if these are present in the initial conditions.] In the text we use the abbreviated 

notation E^ for E V £ m , as we have done with the labelling of the tensor harmonics 

themselves. The summation and integral in QA.27 ) defines the symbolic summation 

over harmonic modes, ^fc; employed in the text. In a closed universe, equation ( |A.27| ) 

should be replaced by 

oo v— 1 £ 

E ab = E \ K \ 3/2 » 2 E E (EvimQ^ + E uim QM). (A.28) 



v=3 



1=2 m= 
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Under the assumptions of statistical homogeneity and isotropy, the covariance structure 
of the E k and E^ is restricted to the form 

(E k E k /) = E {v)b~kk> 



(E k E k ,) = E 2 (v)5, 



kk> 



(E k E k/ ) = 0, (A.29) 

where the symbolic S kk ' represents 5uiS mm i\K\^ 3 / 2 i i ^ 2 S(u — v') in an open universe. In 
the closed case 5{y — v') should be replaced by 5 UU > in this expression. 
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